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SUMMARY 

Samp led- data analysis and harmonic balance analysis are applied to analyze switching DC-DC 
converters under constant on-time control. Design-oriented boundary conditions for the period- 
doubling bifurcation and the saddle-node bifurcation are derived. The required ramp slope to avoid 
the bifurcations and the assigned pole locations associated with the ramp are also derived. The 
derived boundary conditions are more general and accurate than those recently obtained. Those 
recently obtained boundary conditions become special cases under the general modeling approach 
presented in this paper. Different analyses give different perspectives on the system dynamics and 
complement each other. Under the sampled-data analysis, the boundary conditions are expressed 
in terms of signal slopes and the ramp slope. Under the harmonic balance analysis, the boundary 
conditions are expressed in terms of signal harmonics. The derived boundary conditions are useful 
for a designer to design a converter to avoid the occurrence of the period-doubling bifurcation and 
the saddle-node bifurcation. 
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1 Introduction 



Many efforts liave been made in tlie past three decades to model the switching DC-DC converter 
under fixed frequency control [H [H [3l [H [5] . Fewer efforts [H [71 [8l [9] have been made to model the 
switching DC-DC converter under variable frequency control. Constant on-time control (COTC) is 
a type of variable frequency control. In two recent references [SIE], converters under COTC, with 
the inductor current and the output voltage respectively as the feedback variable, are analyzed 
separately. In this paper, these two types of control schemes are analyzed in a unified model. 
The methodology proposed in this paper is based on the sampled-data analysis or the harmonic 
balance analysis, which provides an alternative way besides the tradition approaches to analyze 
the DC-DC converter under COTC. In [HI |9], approximate analysis is directly applied to analyze 
the converter. In this paper, exact sampled-data analysis or harmonic balance analysis is applied. 
Both the sampled-data analysis and the harmonic balance analysis produce the same results, but 
give different and complementary perspectives about the converter dynamics. 

Two bifurcations commonly seen in DC-DC converters are analyzed. They are period-doubling 
bifurcation (PDB) and saddle- node bifurcation (SNB). The boundary conditions associated with 
the bifurcations are derived. Note that, here, the boundary condition means the critical condition in 
the converter parameter space, not about the well known critical eigenvalues in the complex plane. 
The boundary conditions define the bifurcation boundaries in the parameter space to separate 
stable and unstable regions. The methodology proposed here can be applied to general DC-DC 
converters, and the buck converter is used as an example throughout the paper. It will be shown 
that the boundary conditions obtained in [71 El E] become special cases in terms of the general 
models used in this paper. Although the sampled-data analysis has been applied to obtain similar 
boundary conditions in [7], the approach applied here is different. In [7], an approximate sampled- 
data model is directly applied, whereas here an exact sampled-data model is directly applied and 
further approximations are then applied later to preserve the accuracy. 

The remainder of the paper is organized as follows. In Section 2, the operation of COTC is 
reviewed. In Section 3, exact sampled-data analysis is applied. In Sections 4 and 5, PDB and 
SNB are respectively analyzed based on the sampled-data analysis. In Sections 6 and 7, PDB and 
SNB are respectively analyzed based on the harmonic balance analysis. Conclusions are collected 
in Section [51 

2 Operation of Constant On-Time Control (COTC) 

The operation of a switching DC-DC converter under COTC can be described exactly by a unified 
block diagram model |101 [TT] shown in Fig. [H Denote the control signal as Vc which controls the 
inductor current or the output voltage Vq as discussed below. Denote the source voltage as Vg- 
In the model, Ai,A2 G R^^^, Bi,B2 S R^^^^ C,Ei,E2 S R^^^, and D G K^''^ are constant 
matrices, where A'^ is the system dimension. In the n-th cycle, the dynamics is switched between 
two stages. Si and 82- Switching from to S2 occurs at t = J27=i + (where d is fixed in 
COTC). Switching from S2 to 5*1 occurs at t = Tj (where Tj varies in each cycle) when the 

ramp signal h{t) intersects with the signal y := Cx + Du £ R. 

In this paper, two COTC schemes are analyzed. The first scheme is the valley current control 
(a variation of current mode control but with COTC, see Fig. [2] for a buck converter). One has 
y = RiiL — Vc, where Ri is the sensing resistance. At the switching instants, one has y{t) = h{t), 
or equivalently, RiiL = Vc + h{t). The converter changes from S2 to Si when RiiL drops below 
Vc + h{t). The second scheme is the valley voltage control (a variation of control but with 
COTC, see Fig. [3] for a buck converter). One has y = Vq — Vc- The converter changes from ^2 to 
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Figure 1: Block diagram model for switching converter under variable frequency control. 




ramp h{t) 



Figure 2: A buck converter under C-COTC. 



when Vo drops below Vc + h{t). In this paper, these two COTC schemes are denoted as C-COTC 
and V-COTC, respectively. In C-COTC, the inductor current is the feedback signal, and it is also 
called current mode COTC in [8]. In V-COTC, the output voltage is the feedback signal, and it is 
also called COTC in [9]. 

Denote the ramp slope as iria- In the n-th cycle, denote the ramp amplitude as V/i„ = rriaTn. 
The short notation V/i„, instead of V/i^„, is used for brevity. This applies to other variables. Since 
T„, varies for each cycle, the ramp amplitude Vh also varies for each cycle. Let the steady-state of 
the period be T. In steady state, Vhn = Vh = fngT. Denote the switching frequency as fs ■= 1/T 
and let LOg '■= "iirfs- In COTC, the switching occurs at t = Y17=i ^« + where d is fixed, and the 
period T„ is used as a control variable. This is different from the fixed frequency control, where 
the switching occurs at t = nT + dn where dn is controlled, and the switching period T is fixed. 
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Figure 3: A buck converter under V-COTC. 

3 Exact Sampled-Data Dynamics of a General Converter 

3.1 Large-Signal Analysis 

In the fixed frequency control, the states are sampled at t = nT with a fixed period. In COTC, 
the states are sampled at t = Ti with a variable period. Let n„ = [vsmVcn)' ■ Similar to the 

analysis in the large-signal nonlinear constrained dynamics, mapping from x„ := x{^^~^Ti) 
to Xn+i, is 

= e^2(^"-'^)(e^i'^x„+ / e'^^''daBiUn)+ / e^^^^-'^daSs^Xn (1) 

Jo Jd 

g{Xn,Un,Tn) = Cf{Xn, n.„, Tn) + DUn - niaTn 

= (2) 

where the constraint g{xn,Un,Tn) = determines the switching instant T„. In the n-th period, 
once the variable period T„ is determined, Xn is mapped to Xn+i through the function /. 

3.2 Steady-State Analysis 

The periodic solution x^{t) of the system in Fig. [1] corresponds to a fixed point x^{0) in the sampled- 
data dynamics. Let x^{0~) = ^2X^^(0) + B2U denote the time derivative of x^{t) at t = 0~. Let 
y^{t) = Cx^(t) + Du. In steady state, y^{t) = Cx^{t). Let the steady-state duty cycle be D, then 
d = DT. Confusion of notations for the capacitance C and the duty cycle D with the matrices C 
and D can be avoided from the context. 
In steady state, 

pd 

x\d) = e^"^x"iO)+ / e^'^'daBiu (3) 

Jo 

[■T~d 

xO(0) = e"^2(^-'^)x°(d) + / e^2"da52U (4) 

JO 

From ([3|) and dH, one has 

xO(0) = (/-e^2(^-'^)e^i'^)-ne"^^(^-'^) [\'^'''daBiu + \^'''daB2u) (5) 

Jo Jo 
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Let Bi := [Bii, B12] and B2 ■= [B21, -B22] to expand the matrices into two columns. The COTC 
buck converter generahy has Ai = A2, B21 = 0, and B12 = i?22- Then, from ([3]) and ([3]), one has 



^O(o) = (/ _ e^^^)-'A^\e^^^ - e^^(^-'^))Sn^^. - A^'B^2Vc (6) 

Generally the controller may include an integrator (with a pole at zero), making or / — 
^A2(T-d)^Aid non-invertible. In that case, the pole at zero can be replaced by a very small number 
6, then Ai or I — e^2[T-d) ^Aid invertible. Therefore, the invertibility of Ai or I — eA2{T-d)^Axd 
is not critical and can be resolved. This statement about invertibility of a matrix is not repeated 
later. 

The boost converter generally has Bi = B2, then 

xO(0) = (/ - e^2(^-'^)e^^'^)-^(e^^(^-'^) [\^^'^da + r~\^^"da)B^u (7) 

JO Jo 



3.3 Small-Signal Analysis 

Using a hat " to denote small perturbations (e.g., Xn = Xn — x^{Q)). From [21 [lOl [TTl [12], the 
linearized sampled-data dynamics is 



Xn+l = ^Xn + VUn = ^Xn + TiVgn + r2'0cn 

where $ € R,^^^ and V = [ri,r2] € R^^^ (expanded into two columns) are 

df df dg 1 dg 



(8) 



dx 
{I- 



dTn dTn dXn 

x°{0-)C 



{x„,u„,T„)={xO{0),u,T) 



df _ df dg ^1 dg 

dUn dTn^dTj dUr. 



^A2{T-d)^A,c 



(Xn,Un,Tn) = {xO(0),U,T) 



(9) 



(10) 



Local orbital stability of the converter is determined by the eigenvalues of <I>, which are also the 
sampled-data poles of the converter . The periodic solution x^{t) is asymptotically orbitally stable 
if all of the eigenvalues of <1> are inside the unit circle of the complex plane. 

The PDB occurs when one eigenvalue of $ is —1, and det[/ + ^] =0. The SNB occurs when 
one eigenvalue of <^ is 1, and det[7 — <1>] =0. Since other bifurcations [13], such as pitchfork or 
transcritical bifurcations, also have one eigenvalue of ^ at 1, the same analysis can be applied to 
analyze these bifurcations, omitted to save space. All slope-based boundary conditions derived in 
this paper are based on the closed form expression of ([9]). 



3.4 Control-to-Output and Audio-Susceptibility Transfer Functions 

In (jlOp . with the matrices Bi, B2, and D replaced by their first columns, then T becomes Ti. 
Similarly, with the matrices Bi, B2, and D replaced by their second columns, then T becomes r2. 
Since the matrix D = [0, —1] for C-COTC and V-COTC and the second columns of the matrices 
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Bi and B2 are generally zeros, the matrices Ti and r2 can be further simplified as 

r, ^ il - J^^^J ie^^^-^ f e^-^^i^n + f/e^-daB..) (11) 

- CiO(0-)-m, ^^^^ 

Since the output voltage may be discontinuous (as in the boost converter), let E := {Ei + E2) /2 
so that Ex is the average output voltage. From ([8]), the control-to-output transfer function is 

Toc{z) = ^ = E{zl - ^)-^T2 (13) 

The DC gain is Toc{e^^) = 7oc(l) = -E(/ — ^)~^T2- Given a transfer function in the z domain, 
say T{z), its effective frequency response [HI p. 93] is T{e^'^'^), which is valid up to half the 
switching frequency. 

Also from the audio-susceptibility is 

Tosiz) = ^ = E{zl - ^r'Ti (14) 
Other transfer functions, such as output impedance, can be similarly derived as in [TO] . 

3.5 General Slope-Based Boundary Conditions for PDB and SNB 

Theorem 1 Suppose that X is not an eigenvalue of e^'^^'^~'^^ e^^*^ . Then A is an eigenvalue of ^ if 
and only if 

C{I - A-ie^i'^e^2{T-d))-i^O(o-) ^ ^-^^-^ 
Proof: Suppose A is not an eigenvalue of e^^iT-d) ^Aid ^ then 
det[A/-«>] = det[A/-e^2(T-d)gAidj . 

det[I+(A/-e^^(^-'^)e-"^'^)-^( ^.^°^°"^^ )e^2(T-d)^A,d 
= det[A/-e^^(^-'^)e^^'^]- 

(1 + Ce'''(T-d)^Md^xi - e^^(^-'^)e^^'^)-i( ^ )) 

Therefore, A is an eigenvalue of $ if and only if 

1 + Ce^^(T-d)^A,d^^j _ ^MiT-d)^A,dyi( ^°(0-) ) ^ 

Cx^(0~) — nia 

which can be rearranged as 

Cx°(0-) + C(Ae-^2(T-d)g-Aid _ ^^-i^O(Q-) ^ ^^Q^ 

leading to ^ based on the matrix equality /+(Ae~^2(T-d)g-Aid_j-)-i ^ _ x~i^Aid^A2{T~d)yi _ 

□ 
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Note that the condition (|15p is apphcable to general switching converters of any system dimen- 
sion. Also note that in (fT5|) . the left side is related to the ripple slope x^{0~), and the right side is 
the ramp slope nia- As in the popular slope-based boundary condition for the subharmonic oscil- 
lation in the current-mode control, the condition (llSp is also slope-hased. For designation purpose, 
the left side of (jlSp is called an "S plot", S{X,D), as a function of A and D = d/T, for example. 
The S plot can be a function of other variables, such as T, and the power stage parameters, such 
as R, L, and C. 

Corollary 1 

(i) If the system parameters correspond to an occurrence of period- doubling bifurcation (X = —1), 
then 

S{-1, D) = C{I + e^^™e^^^(i-^))-ii;°(0-) = (17) 

(ii) If the system parameters correspond to an occurrence of saddle-node bifurcation (\ = 1), then 

5(1, D) = C{I - e^^™e^^^(i-^))-ii;°(0-) = m„ (18) 

(Hi) Without the ramp (nia = 0), then det[$] = 0, and <1> has an eigenvalue of 0. Equivalently, 
S{0,D) = 0. 

Proof: Proof for (i) or (ii) is directly from Theorem 1. Proof for (Hi) is as follows. 
From Q with = 0, then det[$] = det[/ - x°(0"")C/(Ci;°(0~))] det[e^2(T-d)gAidj ^ ^^^^^ _ 
Cx°{0-)/{Cx°{0-))] det[e^2(T-d)gAidj ^ Another proof for (Hi) is based on the fact that A = 
is a root of (jl6p with rua = 0. □ 
A sampled-data pole at generally causes deadbeat effect [15] where transient perturbation 
correction occurs within N switching periods [?]• A sampled-data pole at also generally causes 
the continuous-time dynamics to have a quality factor Q = 2/7r [3]. 

3.6 Pole Locus and Pole Assignment 

In the traditional root locus plot, the root locus is a function of the feedback gain. In COTC, the 
ramp slope rUa determined by ([15]) can be expressed as a function of the real pole A. Without the 
ramp {nia = 0), one pole is always zero. With the ramp, the poles are shifted according to (jl5p . 
For designation purpose, the S plot as a function of the real pole A alone, is called a pole locus plot. 
Given a desired real pole location, one can see the required ramp slope for pole assignment based 
on the pole locus plot ([15]). 

In summary, the S plot is useful for design purpose to avoid PDB and SNB instabilities, and 
to assign the real poles. First, the S plot S{—1,D) can predict the occurrence of PDB and the 
required ramp slope to avoid PDB. For example, given a ramp slope m-a, the intersection of the 
S plot and the horizontal line at nia shows the unstable operating range of the duty cycle D if 
5(A, D) > rUa- To avoid PDB, one need to increase the ramp slope rUa such that <S'(A, D) < ma for 
the operating range of D. Similarly, the S plot S{1,D) can predict the occurrence of SNB and the 
required ramp slope to avoid SNB. 

Second, given a duty cycle D, a special case of the S plot, S{X,D), as a function of A alone 
(also called pole locus plot), can predict the real pole locations and the required ramp slope for 
pole assignment. For example, the intersection of the S plot and the horizontal line at ma shows 
the real pole location. As the ramp slope varies, the pole location also varies. One can see the 
required ramp slope to assign the real pole to a particular location. If the S plot does not intersect 
with the horizontal line at ma, there are no real poles. Complex poles may exist in that situation. 
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3.7 Approximate Pole Locus Plot 

The exact pole locus plot S{X,D), the left side of (jlSp . can be expressed in an approximate form. 
Based on the assumption that RC and \fLC are much larger than T, matrix approximations such 
as e^^'^ ^ 1 + AiT and (/ + AiTy^ ^ I - AiT can be applied. Then, ^ leads to 

« ^— + )x(0) = ma (19) 

Without the ramp {ma = 0), (fT9]) leads to 

For the buck converter with A\ = A2 and i?2 = 02x2, then i;'^(0~) = Aix''(0) and (f20l) becomes 

c^frx°(o) 

^ - ^ " C^ixO(O) ^^^^ 

where the exaci expression of x''(0) can be obtained from ([5]). Therefore, without the ramp (jria = 
0), two poles of the buck converter are and ([2T]) . Let the state be x = {iiiVc)' , where ii is 
the inductor current and vc is the capacitor voltage. Then, for the buck converter, x'^(O) can be 
approximated as [vgD/R — VsD{l — D)T/2L,VsD]' based on the traditional average analysis jl5j . 

4 Sampled-Data Analysis of PDB in the Buck Converter 

The analysis above is for the general DC-DC converter. The general analysis is applied below to 
the buck converter. Many critical parameters (such as the ramp slope and the constant on-time) 
at the stability/instability boundary can be derived based on the boundary conditions. 

Using the boundary condition to determine the minimum ramp slope. 

The buck converter generally has Ai = A2, B21 = 0, and B12 = i?22- Using ([6]), the PDB boundary 
condition (fTTl) becomes 



Si-1,D) = C{I - e^^-^)-\e^-^ - e^^^(^-^))i?n^. = (22) 

Generally, with larger ma or smaller v^^ the converter is stable. Therefore, the converter is stable if 

S{-\,D)<ma (23) 

Based on the assumption that RC and \J LC are much larger than T, matrix approximations 
such as e^i'^ « I ^ AxT ^ MT'^ and (I + ^iT)"! I - A^T can be applied. Then, the boundary 
condition (p3]) leads to 

5(-l, D) « i-^CBu + {^^)CAxBuT + ( ^ V ^i^n^')^. < rua (24) 

For CAiBn > CAjBuT (or equivalently, RC > T, in V-COTC, for example), ([MD leads to 

D 

{--CBu + {—)CAiBnT)vs < ma (25) 
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or equivalently, using the fact that d = DT, 



2 ^2 



CBii) < ma 



(26) 



which shows the required ramp slope ma to stabihze the converter to avoid PDB. The required 
ramp slope is linearly proportional to D. Therefore, to operate at a higher duty cycle, a larger 
ramp slope is required. Also note that Vq ~ VsD for the buck converter. Based on (I26p and given 
a fixed Vo, the required ramp slope to avoid subharmonic oscillation is the same, independent of 
Vs- This indicates that COTC has very good line regulation to avoid the subharmonic oscillation. 
However, note that the condition (j26p is approximate under the assumption that T is small, the 
exact condition is (|22p . In COTC, d is fixed and T may vary to a large degree. For a large steady- 
state value of T, ()24p is more accurate than ()25p . However, (I24p is complex. One can use ()25p for 
simplicity at the expense of accuracy. 



Using the boundary condition to determine the maximum constant on-time. 

The boundary condition (|26p can be expressed in terms of the maximum constant on-time d, 



d CBn+^ 
2 CAiBii 



if CAiBii > 



(27) 



The inequality sign is reversed if CAiBn < 0. 

Added with external ramps, the two control schemes, V-COTC and C-COTC, with or without 
external ramps added, are analyzed as follows. For V-COTC, one can also use the inductor current 
as the ramp instead of the external ramp h{t). 



4.1 V-COTC Added With an External Ramp (m„ 7^ 0) 



In V-COTC, y 



Vo — Vc- let the load be R, the inductance be L, the capacitance be C, and the 

Let p = R/{R + Rc). For ij^ = 0, p = 1. Then, 



equivalent series resistance (ESR) be R, 

Al : 



Ao 



L 
1_ 

C 



^1 
L 

-1 

RC 



B 
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El = 

p[Rc, 1] 



C = 
CBn 

The boundary conditions (p6]) and ([2 



L 





E2=p[Rc 1 
D=[0,-1] 



pRc 
L 



_2i 
LC 



^(1 



m„ 



d 
2 



> 



< 



CAiBn 
become, respectively, 

R^C s RnC 



L ' 



Dvsf^ d 
2LC 4^ 

pRc I 2ma 
L Dv. 



L 



P 



^(1 



R^. 
L I 



which are further simplified to (if R^C <^ L and R^ <^ R) 



rria 

d 
2 



> 



Dvs ,d 
2LC^2 



< RcC + 



RcC) 
2maLC 



(28) 

(29) 
(30) 

(31) 

(32) 
(33) 

(34) 
(35) 
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Generally the required ramp slope condition is normalized by sj := RcVgD/L, the off-time 
inductor current slope multiplied by Rc- Note that s/ is also the output voltage ripple slope 
contributed by the inductor current. The required ramp slope (|32p becomes 

ma ^ P% d RlC 1 

Sf Z ZRcC L p 

For RlC < L and Rc < R, ([MD is further simplified to 

which agrees with O Eq. 21]. Note that ([37|) is a special case of the general condition pT|) 
(applicable to all kinds of converters). Also note that if the conditions R^C <^ L and Rc ^ R are 
not met, ([55]) will be shown to be more accurate than the boundary condition ([57]) as obtained in 

4.2 V-COTC Without an External Ramp (m„ = 0) 

Without the ramp (ma = 0), the boundary condition (j33|) becomes 

d RpC , , 

2 p(l-¥) 

For ^ ^ and Rc < -R, ([38]) is further simplified to 

^ < RcC (39) 

which agrees with [Hin]. Note that if the conditions RcC ^ L and Rc <^ R are not met, ([38]) will 
be more accurate than the boundary condition ([39|) as obtained in [TJ [9]- 

Pole Locus 

The pole location can be also expressed in a closed form. With some approximations for RcT <^ L, 
(|2T]) leads to 

By setting A > — 1 in (j40p , the maximum of on-time to avoid subharmonic oscillation is 

d RcC^^^j^^^ 



Based on some examples, (j^T|) is more accurate than (j38|) and p9|) . Note that (jlT|) depends on T, 
whereas ([Ml) and are independent of T. For T <^{R + Rc)C, dH]) becomes ([39]) . 
For T - d <C and i?c < -R, (jlQ]) is further simplified to 

2T 
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Figure 4: The pole locus plot ^(A, 0.4) shows the required ramp slope as a function of the pole A. 
Its intersection with the horizontal line at = shows the poles at and -1.1, (solid line for the 
exact expression (jlSp and dashed line for the approximate expression (|19p ). 



which agrees with [3 p. 2679]. Note that if the conditions T — d <^ RC and Rc <^ R are not met, 
(Uni) will be more accurate than the pole location ()12|) as obtained in [7]. This sampled-data pole 
can be mapped to an equivalent continuous-time pole, 

(43) 



2R,C + T-d R^c+^^ 
which is close to the pole l/RcC derived in [H [16] if RcC > T. 

Example 1. [Without compensating ramp.) Consider a V-COTC buck converter from [7] with 
the following parameters: = 5 V, T = 3 //s, d = 1.2 //s, R = 0.5 0, Rc = 20 voTl, L = 2 /iH, 
and C = 20 Here, D = d/T = 0.4. With these parameters, subharmonic oscillation occurs as 
shown in [TJ Fig. 9]. 

Based on the exact sampled-data analysis by calculating the eigenvalues of one pole is 
and the second pole is -1.1. Based on ()40p . the second pole is exactly -1.1. For comparison, based 
on (j42p as derived in [7], the second pole is -1.3. The error is due to the fact that the condition 
T — d <^ RC is not met. The exact pole locus plot 5(A, 0.4) based on (|15p and the approximate 
one (I19p are shown in Fig. [H The intersections of the pole locus plot and the ramp slope (m^ = 0) 
are the two poles, -1.1 and 0. The approximate pole locus plot matches very well with the exact 
one, both accurately predicting the two poles. 

The pole locus plot shows two types of information. First, it shows the resulting pole locations 
given a ramp slope ma- For example, draw a horizontal line at ma in the pole locus plot, the 
intersections with the curve (fTSj) (or (fT9l) ) are the pole locations. Second, given a desired pole 
location, the pole locus plot shows the required ramp slope to assign the pole to that particular 
location. For example. Fig. [5]shows that if the ramp slope is around 9500, the poles can be assigned 
at -0.5 and -0.2. □ 

Example 2. (Pole assignment.) Continued from Example [H let ma = 9500. Based on the exact 
sampled-data analysis by calculating the eigenvalues of the poles are exactly -0.5 and -0.2 as 
predicted. □ 

Example 3. (Determine the minimum ramp slope for stabilization.) The ramp slope does 
not need to be as large as 9500 for stabilization (since only the poles being within the unit circle 
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Figure 5: The S plot «S'(— 1, D) shows the required ramp slope to avoid PDB. The S plot also shows 
that, without a ramp {ma = 0), PDB occurs aX D = 0.36, (solid line for the exact S plot ([22]) and 
dashed line for the approximate S plot (IH 



are required). Based on the exact sampled-data analysis by calculating the eigenvalues of <1>, the 
minimum ramp slope for stabilization is 943.4. From ()22p . the minimum ramp slope for stabilization 
is also exactly 943.4. □ 

Example 4. {Determine the minimum ramp slope for a range of source voltage to maintain 
stability for line regulation.) In Example [TJ the source voltage is fixed at 5 V. Now suppose the 
converter is designed to operate in a range of source voltage Vg G [2, 10]. With a fixed d = 1.2 ^s 
and a fixed Vo = 2, then D G [0.2, 1], T G [1.2, 6] /is, and Us G [1.047, 5.236] rad/s. Based on the 
exact sampled-data analysis, the converter is unstable for D > 0.36 (In Example [H D = 0.4 > 0.36 
and the converter is unstable). 

The required ramp slope based on the exact S plot S{—1,D) (p2|) and the approximate one 
based on (|24p are shown in Fig. [5l The approximate one matches well with the exact one. The 
S plot shows exactly that the converter is unstable for D > 0.36 if no ramp is applied (m^ = 0). 
The S plot also shows that a ramp slope greater than 4217 is required to maintain stability for this 
operating range. □ 

Example 5. {Determine the PDB point.) Continued from Example [T] without a ramp, based on 
Fig. [5l PDB occurs at -D = 0.36. With a fixed d and a fixed Vq = 2, then PDB occurs at T = 3.33 
/is and Vs = 5.56 V. With these parameters and based on the exact sampled-data analysis, the 
eigenvalues of ^ are exactly and -1 as predicted. □ 

Example 6. {Determine the maximum on-time.) Continued from Example [H by calculating 
the eigenvalues of the maximum on-time d to avoid subharmonic oscillation is 1.06 f-ts (with 
T = 2.65 /is to keep the duty cycle at 0.4). The maximum on-time based on (|38|) . (|39|) and (|4ip 
are 0.84 fis, 0.8 fis, and 1.077 /US, respectively, and ()4ip gives the closest prediction. □ 



4.3 V-COTC Added With the Inductor Current as a Compensating Ramp 

Here, the ramp h{t) is not used, and = 0. Instead, the inductor current is sensed through a 
resistance Ri and added with the output voltage to compared with Vc to determine the duty cycle 
as shown in Fig. [6l 
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Figure 6: A buck converter under V-COTC with the inductor current as a compensating ramp. 

The converter is turned on for a fixed time duration d. The converter changes from the OFF 
stage to the ON stage when Vq + iiRi drops below Vc- One has y = Vo + ihRi — Vc and 

C = p[Rc^] + [R^M i^ = [0,-l] 
The boundary condition ([27|) becomes 

d. {pRc + R^)C 

P \^ L pL ) 

For Rc{Rc + Ri)C <^ L and Rc < -R, ([Ml) is further simphfied to 

^<{Rc + Ri)C (45) 

which agrees with [3 Eq. 16]. Note that if the conditions Rc{Rc + Ri)C <^ L and Rc <^ R are not 
met, (|44p will be more accurate than the boundary condition (|45p as obtained in [9], and (j45p is a 
special case of (HH). 

The pole location can be also expressed in a closed form. With some approximations for 
RcT < L, dm leads to 

, 2pr T-d R^{T-d) 

2{Rc + Ri)C + T-d^2RC 2L ' ^ ' 



By setting A > — 1 in (j46p . the minimum of i?j to avoid the subharmonic oscillation is 

2d - iRcC - ^^^^^ , , 

> 4C+^ '''' 

Based on some examples, (j47p is more accurate than (I44p and (j45p . 

Example 7. [Better 'prediction of minimum Ri.) Continued from Example [H by calculating the 
eigenvalues of ^, the minimum Ri to avoid the subharmonic oscillation is 1.82 mi}. The minimum 
Ri based on (jl5|) and ([T7|) are 8.4 mfi, 10 mi7, and 3.4 mO, respectively, and ([T7|) gives the 
closest prediction. □ 



14 



4.4 C-COTC 



In C-COTC, one has y = Riii — Vc and 

C=[Ri,0] D = [0,-1] 



CBii = ^ CAiBn = ^^^^ 



The boundary condition (|26p to avoid SNB becomes 

DvsRi ,dpRc 
2ir^ 2L 



+ l)<ma (49) 



The inequality is always satisfied even without the ramp [rua = 0), and PDB does not occur, agreed 
with IE]. 

Also based on the pole location, one can prove that PDB does not occur in C-COTC as follows. 
With some approximations for R^T <^ L, (j2ip leads to 

Generally <^ LC, then A ~ 1 and PDB (which requires A = — 1) does not occur. This sampled- 
data pole can be mapped to an equivalent continuous-time pole, 

which will be shown to be more accurate than the pole [2L/R — d)/2LC as derived in [8]. 

Example 8. [Better prediction of the continuous-time pole.) Consider a C-COTC buck converter 
from [8] with the following parameters: Vg = 13.2 V, Vo = 3.3 V, T = 1.04 ^us, d = 0.26 /is, i? = 10 
J), Rc = 4.5 mO, Ri = 150 mTL, L = 3.1 jM, and C = 300 /xF. Here, D = d/T = 0.25. 

Based on the exact sampled-data analysis by calculating the eigenvalues of one pole is 
and the second pole is 0.9995 (corresponding to a continuous-time pole at 473 Hz). Based on 
([5T]) . the second continuous-time pole is 419 Hz. For comparison, based on [8], the second pole 
is {2L/R — d)/2LC = 194 Hz. This example shows that (j5ip gives a better prediction of the 
continuous-time pole. 

The exact pole locus plot ^(A, 0.25) based on (fTSj) and the approximate one based on (fT9ll are 
shown in Fig. [71 The approximate one matches very well with the exact one (except that the exact 
one has a finite ramp slope), both accurately predicting the two poles around and 1. The pole 
locus plot also shows that if the ramp slope is small and negative, the two poles may be unstable 
(with the first pole smaller than -1 and the second pole greater than 1), which will be discussed in 
Example □ 



5 Sampled-Data Analysis of SNB in the Buck Converter 

Occurrence of SNB in COTC has not been reported. As discussed above, for V-COTC with a 
negative ramp slope, SNB may occur. In steady state, 

y°(r) - h{T) = Cx°{T) + Du- h{T) = (52) 

which is an equation in terms of T. Generally, if this equation has multiple solutions of T, SNB 
may occur if a converter parameter varies. When SNB occurs, only a single solution exists, which 
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Figure 7: The pole locus plot S{X, 0.25) shows the poles around and 1 without the ramp, (solid 
line for the exact expression (jlSp and dashed line for the approximate expression ()19p ). 



means that the curve Cx^iT) + Du — h{T) has a flat slope (or equivalently, the curve Cx^iT) + Du 
is tangential to h[T)) when SNB occurs. An SNB is also called a tangent bifurcation. 

Take the buck converter for example. The buck converter generally has Ai = A2, B21 = 0, 
and B12 = i?22- Using ([6|), and differentiate ([52]) with respective to T, one has the SNB boundary 
condition, 

5(1, D) = C{I - e^^^)-2(e^^^ - e^''^^'-''^)BnVs = (53) 

One can derive or prove ([53]) in another way. Using ([6|) , the boundary condition (fT8|) directly leads 
to (|53l) . It is interesting to know that the boundary conditions derived based on the steady-state 
analysis and based on the stability (eigenvalue) analysis are the same. 

Based on the assumption that RC and \JLC are much larger than T, matrix approximations 
such as e^i"^ I + AiT + AIt"^ /2 and {I + AiT)'^ ^ I - AiT can be applied. Then, the boundary 
condition ([53]) leads to 

D 2D^ — D 
5(1, D) « {-CA^'Bn - ^^^n + ( )CAiBuT)vs = (54) 

5.1 V-COTC: No SNB With a Positive Ramp Slope 

Using (|31|) which results in CA'^^Bn = — 1 for V-COTC, the boundary condition (|54p becomes 

D 2D^-DTp\ RlC ma .... 

Generally, RlC < L and < 12LC, and the left side of ([55]) is negative, whereas the right side 
is positive with a positive ramp slope. The boundary condition (j55p cannot be met, and SNB does 
not occur in this situation. 

5.2 C-COTC: No SNB With a Positive Ramp Slope 

Using (|48p which results in CA^'^Bn = —Ri/R, the boundary condition (|54p becomes 

D 2D^-D pR,T rUg 

TR 2L ^ 12 ^ L2 ^ vMi ^ ' 
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Generally, T'^RRc < 12L2, and the left side of ([MD is negative, whereas the right side is positive 
with a positive ramp slope. The boundary condition (|56p cannot be met, and SNB does not occur 
in this situation. 

As discussed above, in C-COTC with = 0, one pole is and the second pole is very close 
to 1. The second pole may be shifted to be greater than 1, then SNB occurs. With a very small 
negative ramp slope, SNB may occur as shown in the next example. 

Example 9. {Pole assignment with a negative ramp slope.) Continued from Example [HI based 
on (|56p . a negative ramp slope < —67445 would result in an unstable pole greater than 1. Let 
ma = —100000, which is — 0.42s/, where Sf := VgDRi/L is the off-time inductor current slope 
multiplied by Ri. Compared with Sf, this ramp slope is small. However, this ramp slope will result 
in SNB as shown below. Based on the exact sampled-data analysis by calculating the eigenvalues 
of the poles are -1.675 and 1.0002 as predicted, also agreed with the pole locus plot shown in 
Fig. [71 Although C-COTC is generally believed to be stable even without the ramp, this example 
shows that a perturbation on the ramp slope to make it negative may make the converter unstable 
with occurrence of SNB. □ 



6 Harmonic Balance Analysis of PDB in the Buck Converter 

Consider a buck converter power stage with a control-to-inductor-current (D-to-ii) transfer func- 
tion Gid{s) and a control-to-output-voltage {D-io-Vo) transfer function G^d{s)- In the converter, 
there is an ON switch and an OFF switch (sometimes substituted by a diode). Let the voltage 
across the OFF switch (or the diode) be (as shown in Fig. [2l for example). In the T-periodic 
mode, the waveform of Vd{t) is a square wave with the high voltage at Vg and the low voltage at 0, 
which can be represented by Fourier series (harmonics), 

oo 

Vd{t)= V c„e^""»* where c„ = ^(1 - e"^"'^^'^) (57) 

n=—oo 

Similarly, in the 2T-periodic mode (when PDB occurs), let two consecutive cycle periods heT — 6 
and T + 6. Then, the switchings (from Si to 5*2 and from 5*2 to Si) occurring within the two cycles 
are at t = d, T — 5, T — 5 + d, and 2T. The 2T-periodic signal Vd{t) can be also represented by 
Fourier series, 

oo 

v,{t) = ^r.e'^ -here c„ = (^)(e-^ - 1)(1 + (-l)"e^) (58) 

n=—oo 

In V-COTC, for example, Vo is the feedback signal from the power stage (whereas in C-COTC, 
iL 'is the feedback signal from the power stage). 

In the converter, some parts are linear (from Vd to y) and some are nonlinear (from y to f^). 
Let the Vd-to-Vo transfer function be G^(s). One has [121 p. 470] 

^' Vs LCil + ^)s^ + {^ + RcC)s + l ^ ' 

Let the compensator transfer function (from Vq to —y (negative sign due to negative feedback)) be 
Gc{s). Let the total transfer function from to —y be G{s). Then, one has G{s) = Gc{s)Gy{s) = 
Gc{s)Gvdis)/vs- The gain G{s) is proportional to the loop gain (based on the average analysis) 

Tis) = 9A^)£lM (60) 



17 



by 

G{s) = ^T{s) (61) 

In the s-domain, y = Vc + Gc(s)(fc - Vo) = Vc + Gc{s){vc - Gy{s)vd) = (1 + GciO))vc - G{s)vd. 
Let Re denote taking the real part of a complex number. Then, from (|57p. the T-periodic solution 
y^{t) (at the output of the compensator) is 

oo 

yO(t) = (i + G,{0)>c- J2 Cne^'^^'Gijuus) (62) 

n=— oo 

oo 

= (1 + G,{0))vc - VsDG{0) - 2Re c,eJ""=*G(j?ia;,) (63) 

n=l 

The intersection of h(t) with the T-periodic solution y^{t) = Gx^{t) + Du determines the duty 
cycle and hence the waveform of Vd{t). By "balancing" the equation y^{t) = h{t) (written in Fourier 
series form) at the switching instants, conditions for existence of periodic solutions and PDB can 
be derived. 

6.1 PDB Boundary Conditions Based on Harmonic Balance Analysis 

In the 2r-periodic mode, one has two switching conditions 

y\T-5) = h{T-5) (64) 
y\2T) = h{2T) (65) 

The basic idea to derive the PDB point is as follows. At the PDB point, a period-one mode 
and a period-two mode coalesce, and one has 6 = 0. Subtracting ([M|) from and setting 6 = 0, 
one can obtain the PDB boundary condition 

oo 

^(E (^-"^-i)(-im^)) = -" (66) 



or equivalently, 

oo 

"^.(Re[E(e-"^ - 1)(-1)"G(^)]) = (67) 

n=l ^ 



Since ([66]) . ([67|) (based on the harmonic balance analysis), and ([22]) (based on the sampled-data 
analysis) are exact PDB conditions, they are all equivalent to each other. One has 

oo 

Si-1,D) = |(ReE(e-^ - 1)(-1)"G(^)]) (68) 

n=l 

which is another expression of the S plot in terms of signal harmonics. 

Generally, Gy{s), Gc{s) and thus G{s) = Gc{s)Gy{s) are low-pass filters. The summation in 
(j67|) can be approximated by the term that involves G{s) with the smallest argument, n = 1 for 
example. Therefore, the S plot has an approximate form, 

5(-l,I))«|(Re[l-(e--^)G(^)]) (69) 
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Example 10. (Both harmonic balance analysis and sampled-data analysis produce the same plot.) 
Continued from Example [H with V-COTC, y = Vo — Vc- Then Gc{s) = —1 and G{s) = — 
The plot of (j66p in terms of nia is exactly the solid line in Fig. [5l indicating that (j66p is equivalent 
to (p2]) (based on the sampled-data analysis). □ 



6.2 Harmonic Balance Analysis of PDB in Terms of the Loop Gain 

Since the gain G{s) is proportional to the loop gain T(s) = G{s)vs/Vh = G{s)vs/maT, ([55]) directly 
leads to the following theorem. 

Theorem 2 Given a closed-loop COTC buck converter with a loop gain T{s) as defined in [6U\). 
PDB occurs when 

oo 

E (e-^-l)(-irT(^) = 2 (70) 

?i=— oo 

The condition (j70p can be expressed in various forms, for example, 

oo 

ReE(e-^-l)(-l)"r(^)] = l (71) 

n=l 

The summation in (I7ip can be also approximated by the term that involves T(s) with the 
smallest argument. Therefore, (j7ip becomes 

Re[(l-e-^)r(^)]«l (72) 

It should be noted that (j72p is only an approximate condition, and the exact condition is (|70p . 
The "LI plot" in the real domain. 

Note that the boundary condition (|7ip is a function of d, Ug, and the loop gain T(s), where T{s) 
is further a function of Vg, rua, the power stage and compensator parameters. Define an LI plot, 
which is a real function, as 

oo 

L,(^.):= E (e-^-l)(-irT(^) (73) 

n=— oo 

Then, SNB occurs when 

Li{lOs) = 2 (74) 

Note that uig = 2tt/T = 2TTD/d, which is linearly proportional to D given a fixed value of d in 
COTC, Li[ujs) can be represented as a function of D instead of ojs- 

The "L2 plot" in the real domain. 

Since in some situations, rUa = when no ramp compensation is used, then T{s) = G{s)vs/Vh 
becomes infinite. In such a situation, define an L2 plot, which is also a real function, as 

oo 

L,(..):= E (e-^-l)(-l)"G(^) (75) 

n=— oo 
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Then, from (1661). SNB occurs when 



L2{u^s) = ^ (76) 

Example 11. {Fig. is also an L2 plot) Continued from Example HI since cjg = 27r/T = InD/d 
is linearly proportional to D, Fig. [5] is also an L2 plot with the horizontal axis scaled by LOg = 
2tt/T = 2TTD/d. □ 

The "H plot": a Nyquist-like plot in the complex plane. 

Let 

oo 

(e-^_l)(_l)«G(^) (77) 

n=l 

Then, PDB occurs when 

ITT) 

^e[H{u,)] = ^ (78) 

For designation purpose, H{iOs) is called an H plot because it is similar to the Nyquist plot in the 
complex plane. Given a desired range of ojs = 2'k/T = 2TTD/d, one can plot H{oJs) according to 
(j77p to determine whether PDB occurs in this range of ujs- 

6.3 V-COTC 

As shown below, the harmonic balance analysis also leads to the same conditions derived by the 
sampled-data analysis. In V-COTC, Gc{s) = -1, and G{s) = Gc{s)Gi,{s) = -G.y{s). Using (f59]) 
and the facts that uJsd/2 = ttD and, for < D < 1, 



^ (l-cos(vrfcD))(-l)^ _ TT^D^ 

p - ^ (79) 

fc=i 

sin(7r/fcL>)(-l)'= ttD 



fc2 



fc=l 

the boundary condition (|67|) becomes 



^K^-^cC)<!^^ (81) 



2LC"2 ' Vs 

which is exactly (j35p . but expressed in a different form. It is interesting to note that both the 
sampled-data analysis and the harmonic balance analysis lead to the same boundary condition, 
providing convincing evidence about the accuracy of the derived boundary condition. 

6.4 C-COTC 

In C-COTC, the feedback signal from the power stage is i^. Similar to ([59]) . the Vd-to-ii transfer 
function is [121 p. 470] 



Gidis) _ (1 + f )Cg + ^ 

LC{1 + ^)s^ + + RcG)s + 1 



G,is) := = . ^;:.^/r;"^^. . . (82) 



Since no extra compensator (except the compensating ramp h{t)) is added in the current loop, 
Gc{s) = —1, and the gain G{s) = Gc{s)Gi{s) = —Gi{s). Since the expression of G{s) is derived, 
the rest of analysis is similar to the case for V-COTC. 
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7 Harmonic Balance Analysis of SNB in the Buck Converter 



Using ([57|) and ([62]). in steady state, 



oo 



j2mr 



y°(r) = (i + Ge(o)K- 



j2n7r 




T 



(83) 



n=— oo 



As discussed above, the curve y^{T) — h{T) has a flat slope (or equivalently, the curve ip{T) is 
tangential to h{T)) when SNB occurs. Differentiating (|83|) with respective to T, one has the SNB 
boundary condition, 



8 Conclusion 

Design-oriented bifurcation boundary conditions are derived for general switching DC-DC con- 
verters under COTC. The derived boundary conditions define the instability boundaries in the 
converter parameter space, and are therefore useful for a designer to design a converter to avoid 
the occurrence of the bifurcations or instabilities. Two bifurcations commonly seen in DC-DC con- 
verters are analyzed. They are PDB and SNB. Two COTC schemes, V-COTC and C-COTC, are 
also are analyzed based on a unified model. The derived boundary conditions are more general and 
accurate than those recently obtained in [3 [HI [9]. The boundary conditions recently obtained in 
[71 El [9] become special cases under the general modeling approach presented in this paper. Since 
the analysis is based on the general block diagram model shown in Fig. [H once a converter (of any 
system dimension under various control schemes) is expressed in terms of the block diagram model, 
the boundary conditions can be readily obtained. 

Both exact sampled-data analysis and harmonic balance analysis are applied. Different analyses 
give different perspectives on the system dynamics and complement each other. Under the sampled- 
data analysis, the boundary conditions are expressed in terms of signal slopes and the ramp slope. 
Under the harmonic balance analysis, the boundary conditions are expressed in terms of signal 
harmonics. Although occurrence of SNB in the converter under COTC is rare, it is shown that 
a small perturbation on the ramp to make the ramp slope negative may cause the occurrence of 



Many new design-oriented plots are proposed. The S plot predicts PDB and SNB instabilities 
and the required ramp slope to avoid them. A pole locus plot, a special case of the S plot, predicts 
the real pole location and shows the required ramp slope for pole assignment. The S plot can be 
expressed in terms of matrices or signal harmonics. The S plot in terms of signal harmonics also 
leads to other design-oriented plots, such as the LI and L2 plots in the real domain, and the H plot 
in the complex plane, to facilitate stability analysis and design. 
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